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ABSTRACT 



Context. We present the first dynamical study of the globular cluster system of NGC 4636. It is the southernmost giant elliptical galaxy of the 
Virgo cluster, and is claimed to be extremely dark matter dominated, according to X-ray observations. 

Aims. Globular clusters are used as dynamical tracers to investigate, by stellar dynamical means, the dark matter content of this galaxy. 
Methods. Several hundred medium resolution spectra were acquired at the VLT with FORS 2/MXU. We obtained velocities for 174 globular 
clusters in the radial range 0.'90 < R < 15.'5, or 0.5 - 9R e in units of effective radius. Assuming a distance of 15 Mpc, the clusters are found 
at projected galactocentric distances in the range 4 to 70kpc, the overwhelming majority within 30kpc. The measured line-of-sight velocity 
dispersions are compared to Jeans-models. 
(— ~) Results. We find some indication for a rotation of the red (metal-rich) clusters about the minor axis. Out to a radius of 30 kpc, we find a 
v,Q roughly constant projected velocity dispersion for the blue clusters of cr a 200 km s~ . The red clusters are found to have a distinctly different 
f*"*) ■ behavior: at a radius of about 3', the velocity dispersion drops by ~ 50 km s~ to about 170 km s which then remains constant out to a radius 
of 7'. The cause might be the steepening of the number density profile at ~ 3' observed for the red clusters. Using only the blue clusters as 
dynamical tracers, we perform Jeans-analyses for different assumptions of the orbital anisotropy. Enforcing the model dark halos to be of the 
NFW type we determine their structural parameters. Depending on the anisotropy and the adopted M/L-values, we find that the dark matter 
O fraction within one effective radius can vary between 20% and 50% with most a probable range between 20% and 30%. A main source of 
uncertainty is the ambiguity of the velocity dispersion in the outermost bin. A comparison with cosmological N-body simulations reveals no 
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J2 ' striking disagreement. 

. , Conclusions. Although the dark halo mass still cannot be strongly constrained, NGC 4636 does not seem to be extremely dark matter domi- 

^* . nated. The derived circular velocities are also consistent with Modified Newtonian Dynamics. 

Key words. Galaxies: elliptical and lenticular, cD - Galaxies: star clusters - Galaxies: halos - Galaxies: kinematics and dynamics - Galaxies: 



* 



ft ' individual NGC 4636 

1. Introduction surements based on the globular cluster luminosity function 

(GCLF) presented by Kissler et al. ( I1994I and Dirsch et 

1.1. Work on NGC 4636 al ( 2 005 ) point towards a larger distance of more than 17 Mpc. 

NGC 4636 is the southernmost bright elliptical in the Virgo Tabled summarizes the basic properties of NGC 4636. 
cluster of galaxies. It is located 10° south of M87, within Although NGC 4636, being an elliptical galaxy, predomi- 

the Virgo Southern Extension. Surface brightness fluctuation nantly consists of an old stellar population, there are some weak 

(SBF) measurements yield a distance of 15 Mpc (Tonry et indications of an intermediate age population: firstly, it is a su- 

al. 2001 1, which we adopt throughout this paper. This value, pernova la host galaxy (SN 1939a: Zwicky 1939 Giclas 1939 1, 

however, has to be considered as a lower limit, since mea- and there are theoretical considerations that this type of super- 

nova originates in intermediate age populations rather than in 

Send offprint requests to: ylva@astro.uni-bonn.de very old ones (Yungelson et al. 1995 and Yoshii et al. 1996 1. 

* Based on observations collected at the European Southern Secondly, far infrared observations by Temi et al. (2003) hint 

Observatory, Cerro Paranal, Chile; ESO program 69.B-0366. at & recent ^ which ^ NGC ^ ^ ^ ^ 

Founded by merging of the Institut firr Astrophysik 
und Extraterrestrische Forschung, the Sternwarte, and the In the 1970s NGC 4636 was the target of several Hi- 

Radioastronomisches Institut der Universitat Bonn. observations. Despite some claims of an extended H i-emission 
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(Knapp et al. 119781 Bottinelli & Gouguenheim ll978> . the deep- 
est observations performed by Krishna Kumar & Thonnard 
(1983l, yielded an upper limit of only 11.5 x 1O 7 M0 to the 
total H i mass of NGC 4636. 

Because of its high X-ray luminosity, NGC 4636 has been 
studied extensively in this frequency range and these observa- 
tions have also been used to infer its dark matter content and 
distribution. 

Forman et al. (1985 1 analyzed Einstein X-ray data from a 
large sample of early type galaxies. For NGC 4636 they found 
extended emission (out to a radius of R max - 44 kpc) and de- 
rived an X-ray luminosity of L x = 1.78 + 0.10 x 10 41 ergs s -1 . 
They quote a mass-to-light ratio M tota \IL^ - 87 for /? max , the 
second largest value in their sample. 

Matsushita et al. ( 1998 1, using ASCA, found very extended X- 
ray emission around NGC 4636 out to 60'. This corresponds to 
a radius of almost 300 kpc in a distance of 15 Mpc. The authors 
suggest that the gravitational halo of NGC 4636 itself termi- 
nates in the region between 10 and 30 kpc. They observe an 
increase in mass beyond 30 kpc forming a halo-in-halo struc- 
ture, whose size is comparable to that of a galaxy group. 
Jones et al. (2002i used the Chandra X-Ray Observatory 
Advanced CCD Imaging Spectrometer (ACIS) to study the X- 
ray halo of NGC 4636 with a spatial resolution of ~ 50 pc. The 
observed X-ray features are so far unique: the images show 
symmetric, arm-like structures which surround a bright central 
region. These structures extend ~ 8 kpc from the galaxy center. 
The leading edges of these arms show changes in brightness by 
a factor of two on scales of just a few arcseconds. While these 
sharp fronts might be suggestive of a merger, the authors favor 
a nuclear outburst as the source of the observed morphology. 
They interpret the arms as the projected edges of paraboloidal 



Table 1. NGC 4636 basic data compiled from the litera- 
ture (l)NED, (2) Tonry et al. (|200T1 . (3) thi s work , (4)R3C, 
(5)Dir sch et al. ( 120051 . (6) Bend er et al. fi994l . (7)Idiart 
et al. d2"0"0"3l . (8) Forman et al. JT9851 . (9) Loewenstein et 
al.d200Tl. 



Other names 



UGC07878, VCC1939 



Position (J2000) 
Galactic coordinates 



121 42™ 50* +02 °4L 17'.' 
/ = 297.75° b = 65.47° 



(1) 



SBF distance modulus 


(m-M) = 30.83 ±0.13 


(2) 


Distance 


D= 15 Mpc 




Scale 


l"=73pc 1'= 4.4 kpc 




GCLF distance modulus 


{m-M) = 31.24 ±0.17 


(5) 


Heliocentric velocity 


fheiio = 906 ± 7kms~' 


(3) 


Hubble type 


E0+ 


(4) 


Ellipticity 


6 = 0.15 


(6) 


Position angle 


PA = 145° 


(5) 


Effective radius 


R, = 101'.' 7 (=7.5 kpc) 


(7) 



Blue magnitude 

Total apparent magnitude 



B T = 10.43 
Tl= 8.70 ± 0.05 



(4) 
(5) 



Metallicity 
X-ray luminosity 
Nuclear X-ray emission 



[Fe/H] = -0.01 dex (7) 

L x = 1.78 ±0.10 xl0 41 ergs/s (8) 
<2.7x 10 38 ergs/s (9) 



shock fronts expanding from the nucleus. Jones et al. speculate 
that the unusually high X-ray luminosity of NGC 4636 might 
be due to the advanced state of the cooling that brought on the 
outburst and may be further enhanced by the shocks. When the 
gas halo returns to hydrostatic equilibrium, the X-ray luminos- 
ity could decline by an order of magnitude or more. 
Loewenstein & Mushotzky (2003 1 used the Chandra X-Ray 
Observatory to measure the hot interstellar medium tempera- 
ture profile of NGC 4636. Assuming hydrostatic equilibrium, 
they determined the total enclosed mass profile for the radial 
range 0.7 < r < 35 kpc. These authors come to the conclusion 
that dark matter constitutes a large fraction (between 50% and 
80%) of the total mass even inside the effective radius. 
Given all these peculiarities, it is of great interest to use globu- 
lar clusters as dynamical tracers in order to evaluate the dark 
matter content of NGC 4636. This work continues our pro- 
gramme on the dynamics of globular cluster systems around 
early-type galaxies. Previous papers on NGC 1399, the cen- 
tral galaxy of the Fornax cluster are Richtler et al. (2004i and 
Dirsch et al. (120031. 



1.2. The NGC 4636 Globular Cluster System 

The globular cluster system of NGC 4636 has been the target 
of three photometric studies. Hanes (119771 used photographic 
plates to assess the number and the luminosity of the GCs. He 
derived a specific frequency (i.e. the number of GCs normal- 
ized to the host galaxy's luminosity (Harris & van den Bergh 
I198U of Sn - 9, which is unusually high for a quite iso- 
lated elliptical galaxy. Most field ellipticals have S N values be- 
tween 3 and 4. The first CCD-study was performed by Kissler 
et al. (1994) who observed a field of 7 '5 x 7 '5 centered on 
NGC 4636 in one band (Cousins V, limiting magnitude of 
24.25). They estimated the total number of GCs to 3600 + 200, 
which corresponds to Sti = 7.5 ± 2.0, hence supporting the un- 
usually high number of GCs. 
The full extent of the GCS together with color information only 



became apparent in the recent study of Dirsch et al. (2005 1. 
They used the MOSAIC camera at the 4-m telescope at the 
Cerro Tololo Interamerican Observatory (CTIO). The bands 
were Kron-Cousins R and Washington C. The final images 
have a field-of-view of 34'7x34'7. The color distribution of the 
GCs shows the familiar bimodal shape. The color of the min- 
imum of the inner sample (C-R=1.55) was chosen to divide 
red (metal-rich) from blue (metal-poor) clusters, frequently in- 
terpreted as two different sub-populations (e.g. Ashman & Zepf 
IJ99"8l Kundu & Whitmore l20TJTl Larsen et al. 20® . For the to- 
tal number of globular clusters Dirsch et al. derived 3700 + 200, 
which agrees well with the previous measurement. The spe- 
cific frequency Sn ranges between 5.8 + 1.2 and 8.9 + 1.2, for 
the GCLF and SBF distance, respectively. An interesting fea- 
ture shows up in the radial distribution of the clusters: While 
the (projected) radial number density distribution of the blue 
clusters can be described by one power-law (n(r) oc r") within 
1'- 8', the density distribution of the red clusters changes the 
exponent within this radial interval. At a radial distance of 7' 
(9' ) for the red (blue) clusters, the power-law exponent drops 
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to about -5 (cf. Table 3 in Dirsch et al. I2005J Such behavior 
has not been observed in any GCS before. This abrupt decrease 
in number density might mark the limit of the GCS and pos- 
sibly that of the galaxy itself. The photometric data by Dirsch 
et al. forms the basis of the candidate selection for our spectro- 
scopic observations described in the following section. 

2. Observations 

The observations were carried out during four nights (May 
7-10 2002) at the European Southern Observatory (ESO) 
Very Large Telescope (VLT) facility at Cerro Paranal, Chile. 
The VLT Unit Telescope 4 was used with the FORS2 (FOcal 
Reducer/low dispersion Spectrograph) instrument equipped 
with the Mask EXchange Unit (MXU). The standard resolu- 
tion collimator used for this program provided a field-of-view 
of 6.'8 x 6'. 8. The detector system consisted of two 4096 x 2048 
red optimized CCDs with a pixel size of 15//m. For our ob- 
servations the chips were read out using a 2 x 2 binning. The 
grism 600B gave a spectral resolution of about 3 A. The spec- 
tral coverage was dependent on the slit position on the mask. In 
most cases, the usable coverage was about 2000 A with limits 
on the red side varying between 5500 and 6500 A. We exposed 
13 spectroscopic masks, the preparation of which is described 
in the next section. Flat fielding was done with internal flat 
lamps. A Hg-Cd-He lamp was used for wavelength calibration. 
The observations are summarized in Table|2] 

2.1. Mask Preparation 

Pre-images of the 13 pointings in 7 fields (see Fig. [Q were 
obtained in service mode in October 2000. Each pointing was 
observed in Cousins V and I filters for 30 seconds. The GC 
candidate selection was based upon the photometric work pre- 
sented in Dirsch et al. (2005 1. Cluster candidates had to fulfill 
the following criteria: they had to be brighter than 22.7 in R and 
the allowed color range was 0.9 < C-R < 2.1. Further, the can- 
didates should exhibit a star-like appearance on the pre-images 
to distinguish them from background galaxies. The ESO FORS 
Instrumental Mask Simulator (FIMS) software was then used 
to select the positions, widths and lengths of the slits. A slit 
width of 1" was chosen according to the sub-arcsecond seeing 
conditions on Paranal. To maximize the number of GC spec- 
tra obtained per mask, we chose to observe our targets and sky 
positions through separate slits of 2" length. Obviously, the 
quality of the sky subtraction then depends on the accuracy of 
the wavelength calibration. As it will be shown, it turned out to 
be satisfactory to obtain radial velocities with the desired accu- 
racy. After positioning the slits for the selected GC candidates, 
the remaining space on the masks (especially in the outer fields) 
was used to include additional objects. Thus, also background 
galaxies and point sources not matching the above-mentioned 
criteria were observed. 

2.2. The Data Set 

The MXU-observations of our thirteen fields have exposure 
times in the range of 3600 to 7200 seconds (see Tabled for a 
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Fig. 1. Positions of the fields overlaid on a 33' x 33' DSS im- 
age centered on NGC 4636. For each mask the entire FORS 2 
6'8 x 6!8 field-of-view is shown. Note that the slit positions are 
confined to the central 43 x 6'8 region of each pointing. 



summary of the observations). To minimize the effects from 
cosmic ray hits, the observation of each mask was divided 
into two exposures - with the exception of Field 1 . 1 for which 
three science images were obtained. In addition to the spec- 
troscopic observations, calibration measurements, i.e. bias, flat 
fields and wavelength calibration exposures were obtained dur- 
ing day time. 

3. Data Reduction 

The following section summarizes the reduction steps per- 
formed to obtain the final sky-subtracted spectra. Most of the 
data reduction was done within the IRAF-environment 1 . 

3.1. Basic Reduction 

The merging of the two CCD-frames was done with fsmosaic 
which is available as a part of the ESO-FIMS software. The 
dimension of the combined images is 2050 x 2076 pixels. For 
cosmic ray removal, the task bclean from the Starlink Figaro- 
package was found to perform best. Biasing and flat-fielding 
was done in the standard way. Since the aim of this study is the 
measurement of radial velocities rather than the measurement 
of line strengths, we did not apply any response function fitting 
or flux calibration. To obtain one-dimensional spectra (aperture 
definition, tracing, and extraction) the IRAF task apall from the 



1 IRAF (Image Reduction and Analysis Facility) is distributed by 
the National Optical Astronomy Observatories, which are operated 
by the Association of Universities for Research in Astronomy, Inc., 
under cooperative agreement with the National Science Foundation. 
I http://iraf.noao.edu/ 1 
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Field 


Center Position 


Exp. Time 


Seeing 


# Slits 


#GCs 


Night 


UT 




(J 2000) 


(sec) 




total 


objects 






(start) 


1.1 


12:42:48.2 


02:39:21.0 


7200 


l'.'l 


111 


79 


39 


2002-05-07 


23:27 


1.2 


12:42:45.8 


02:39:24.2 


4500 


0'.' 9 


98 


61 


47 


2002-05-07 


01:49 


1.3 


12:42:50.8 


02:39:21.3 


5400 


1'.'2 


106 


65 


26 


2002-05-07 


03:17 


2.1 


12:42:17.8 


02:37:27.9 


3600 


l'.'O 


101 


51 


3 


2002-05-07 


05:11 


2.2 


12:42:23.0 


02:37:27.6 


3600 


l'.'O 


99 


48 


2 


2002-05-08 


23:32 


3.1 


12:42:34.2 


02:46:22.4 


3600 


0'.'8 


113 


77 


26 


2002-05-08 


01:21 


3.2 


12:42:28.9 


02:46:20.7 


3600 


l'.'l 


113 


54 


14 


2002-05-08 


02:31 


4.1 


12:42:27.1 


02:53:27.4 


3600 


1'.'3 


98 


44 


2 


2002-05-10 


23:37 


5.1 


12:43:00.6 


02:46:50.3 


3600 


0'.'9 


105 


17 


59 


2002-05-08 


03:43 


6.1 


12:43:15.0 


02:39:33.7 


3600 


1'.'2 


107 


53 


8 


2002-05-10 


01:04 


6.2 


12:43:12.6 


02:39:32.0 


3600 


0'.'9 


110 


55 


16 


2002-05-09 


00:18 


7.1 


12:43:33.9 


02:32:12.1 


3600 


0'.'9 


109 


52 





2002-05-09 


01:31 


7.2 


12:43:29.3 


02:32:09.7 


3600 


0'.'8 


113 


59 


2 


2002-05-09 


02:44 



Table 2. Summary of observations (ESO program ID 69.B-0366(B)). The seeing values were determined from the acquisition 
images taken just before the MXU exposures. 



apextract package was used. Once the aperture positions were 
defined on the science images, they were checked and - where 
necessary - adjusted on the flatfield divided exposures. For the 
2" long slits used for most sources, an aperture size of 4 pixels 
(= 1") was found to yield the best results. 

3.2. Wavelength Calibration and Sky Subtraction 

The one-dimensional arc lamp spectra were calibrated using 
the interactive task identify from the noao.onedspec package. 
The dispersion/wavelength solution was approximated by a 
seventh-order Chebyshev polynomial. The residuals of this fit 
were inspected for outliers, and if necessary the deviant data 
points were excluded from the fit in the next step. The RMS 
errors of the final fits were about 0.04 A. As described above, 
sky spectra were measured through separate slits. In order to 
perform the sky subtraction for a given GC-spectrum, two or 
three adjacent (within 18") of these sky spectra were com- 
bined. This combined spectrum was then subtracted using the 
skytweak task. 

4. Measurement of radial velocities 

The IRAF task fxcor from the noao.rv package performs a 
Fourier cross-correlation on the input object and template spec- 
tra. This program is based upon the technique developed by 
Tonry & Davis (1979 1. For the cross-correlation, a template 
with a high S/N and a spectrum similar to that of a globu- 
lar cluster is required. NGC 4636 itself, however is not suit- 
able due to its high velocity dispersion of about 200kms _1 . 
Therefore, the velocities were measured using a spectrum (S/N 
~ 30) of NGC 1 396 as template. This spectrum was obtained in 
Dec. 2000, during observations of the NGC 1399 GCS, using 
the same FORS2 instrumental setup and has already been used 
as template for investigating the kinematics of the NGC 1399 
GCS (Dirsch et al. l2004l . NGC 1 396 is a dwarf elliptical galaxy 
close to NGC 1399 , and its spectrum resembles that of a metal- 
rich globular cluster. Its velocity dispersion is only 65 ±6km s 



(Wegner et al. 2004). We adopted a heliocentric radial velocity 
of 815 + 8 km s _1 (Dirsch et al. 2004 1 . The velocities obtained 
with this template were checked with a second template spec- 
trum obtained during our observations (object 1.3:56), a bright 
star whose observed velocity was found to be -682 km s with 
respect to the galaxy template. The velocities measured with 
the two different templates agree very well. Since the fxcor un- 
certainty estimates were smaller for the galaxy template, only 
these velocities are used for the subsequent analysis. The ac- 
curacy of the velocity determination depends on the shape of 
the spectrum, and the wavelength range for which the cross- 
correlation is performed. A range of 4700 < A < 5500 A was 
found to be an appropriate choice in most cases. The upper 
limit was chosen to avoid the night sky emission lines found 
around the strong [OI] 5577 A feature. In some spectra how- 
ever, the relatively weak Nitrogen line at 5199 A left substan- 
tial residuals: where this was the case, the wavelength interval 
was adjusted accordingly. For the velocities of the NGC 4636 
GCS , we applied a heliocentric correction of i>hc = - 1 8 km s~ ' ; 
the correction for the NGC 1 396 spectrum is v^ c — - 1 km s~ ' . 
The uncertainty for the velocities is estimated using: 



Ay 



-4 



fxcor template 



where Auf XCOT is the velocity uncertainty output by fxcor, and 
Autempiate is the velocity uncertainty of the template spectrum. 

4.1. Zero Point Shifts 

Since the wavelength calibration exposures for each mask were 
obtained during the day following the observation, the mask 
was moved into the focal plane again. The finite positioning ac- 
curacy of the MXU, and shifts due to instrument flexure under 
gravity - all calibration exposures are obtained with the tele- 
scope pointing at the zenith - result in an offset between the 
science and the calibration images. This offset manifests itself 
as a zero point shift in the observed wavelength of the strong 
telluric [OI] 5577 A emission line. The zero point shifts in our 
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Fig. 2. Histogram of the velocity uncertainties calculated by 
the fxcor-task. The grey histogram shows the values for the 
blue clusters, the black histogram shows the same for the red 
clusters. 



Fig. 3. Velocity uncertainties versus R-magnitude for the 
globular clusters (crosses: blue clusters; open circles: red clus- 
ters). 



data set were found to vary with the position of the slit along 
the y-axis of the mask. To quantify this effect, the wavelength 
of the 5577 A line was measured in all spectra for a given mask. 
A linear function was fit to the data and provided a velocity cor- 
rection for every aperture. These corrections have values in the 
range -25km s -1 < Ai>mxu ^ 25kms~ 1 . 

To measure the systemic velocity of NGC 4636, we used 
the nine spectra of NGC 4636 from mask 1_1, and found a 
heliocentric velocity of 906 + 7kms _1 , In the NASA/IPAC 
Extragalactic Database (NED) 2 one finds seven optical ve- 
locity measurements with uncertainty estimates. The mean is 
917 +31kms~ 1 . Hence, our velocity agrees very well with pre- 
vious measurements. 

4.2. Velocity Uncertainties 

For the GCs, no trend is visible for the metal poor (blue) and the 
metal rich (red) subpopulations: both have very similar error 
distributions (see Fig.|2j. Figure[5]shows the velocity uncertain- 
ties as a function of apparent R-magnitude. As expected, the 
measurements of fainter objects tend to have larger uncertain- 
ties. The uncertainties returned by the fxcor task are computed 
according to the formulae given by Tonry and Davis (1979 1. 
To check whether these uncertainties are a realistic error esti- 
mate, we compared the velocities of those 40 objects for which 
we had two spectra since they were present on two masks. The 
mean difference for the velocity measurements was found to 
be very similar to the mean of the uncertainties given by fxcor. 
Hence, the fxcor uncertainties can be considered a good error 
estimate. The velocities determined for the GCs are listed in 
appendix|X]which is available in electronic form. 

5. Properties of the Globular Cluster Sample 

For every object with a velocity measurement, the color and 
magnitude information was extracted from the photometry- 



2 The NASA/IPAC Extragalactic Database (NED) is operated by 
the let Propulsion Laboratory, California Institute of Technology, un- 
der contract with the National Aeronautics and Space Administration. 



database from Dirsch et al. (2005) via coordinate matching. 
The coordinates, velocities, C - R -colors and apparent R mag- 
nitudes of these objects are listed in appendixlAl For two clus- 
ters, neither color nor magnitude information exists, since they 
happened to be located near a bad row of the MOSAIC-CCD. 

5.1. Defining the Sample 

The first step was to divide our data set into GCs and fore- 
ground stars. A compelling distinction between stars and 
globular clusters would be given if the GCs could be resolved 
on the acquisition images. However, only very few, about three, 
of our clusters are marginally resolved. Guided by Fig.@]which 
shows the distribution of velocity versus color for all measured 
objects, a lower limit of 350kms _I for the radial velocity of 
the GCs was chosen. However, one has to be aware that there 
might be a velocity domain where this identification is doubt- 
ful. Given that the Galactic rotation vector towards the posi- 
tion of NGC 4636 is -78kms , one cannot exclude that fore- 
ground stars have heliocentric velocities of up to 400kms~ 1 . 
Neither can one rule out the presence of GCs with radial veloc- 
ities of more than 500 km s _I relative to NGC 4636. 

The range in velocities occupied by foreground stars and 
globular clusters is -248 < v < 313 kms -1 and 392 < v < 
1441 kms -1 , respectively. The data set comprises 174 globular 
clusters and 175 foreground stars. The basic properties of the 
GC sample are listed in Table[5] The so-defined GCs appear in 
Fig. [5] which shows their distribution around NGC 4636. The 
inner region is well sampled, but due to the sharp drop-off at 
~7' only a handful of clusters is located outside this radius. 

5.2. Color and Luminosity Distribution 

For any further investigation, it is important that the color dis- 
tribution of the GCs with measured radial velocities resem- 
bles that of the entire cluster system. Indeed, the bimodality 
of the color distribution found in the photometric study by 
Dirsch et al. (2005 1 is clearly visible in our much smaller sam- 
ple (see Fig. [6ji. For consistency, the same dividing color of 
C-R = 1.55 was used. Thus, we find 97 blue and 75 red 
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all 




blue 




red 


number of GCs 


174 




97 




75 


mean velocity [km s -1 ] 


904 




911 




887 


min. velocity [kms~'] 


392 




514 




392 


max. velocity [kms -1 ] 


1441 




1441 




1273 


mean vel. uncertainty [km s~ ] 


31 




33 




28 


velocity dispersion [<x(R)kms _1 ] 


203 + 


11 


202 ± 


15 


199 ± 17 


color range C-R 






0.85- 


1.55 


1.55-2.15 


range in projected radius 


0:91- 


15:43 


0:91 - 


15:43 


i:00- li:90 



Table 3. Basic properties of the globular cluster sample 
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Fig. 4. Color vs. radial velocity. All objects with heliocentric 
velocities higher than 350kms _I (dotted line) are regarded as 
globular clusters. The objects with velocities falling below this 
limit are foreground stars. The solid line marks the systemic 
velocity of NGC 4636. 
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globular clusters. The color magnitude diagram (CMD) for the 
GCs is shown in Fig. [6] Although the selection of GC candi- 
dates was obviously biased towards bright clusters, no clusters 
brighter than R — 19 were found. A few clusters are fainter 
than R - 22, but their velocities have large errors (cf. Fig. [3} 
and should be considered with caution. Figure shows the 
histogram of the apparent R-magnitudes for the globular clus- 
ters. It illustrates that our spectroscopic study only probes the 
very bright part of the GCLF, well above the turn-over mag- 
nitude at 23.33 (Dirsch et al. 2005 1. Furthermore, it shows 
that there is no significant difference between the distributions 
of blue and red clusters. The magnitude range of the GCs is 
19.11 > inn > 22.73 in Kron R. When converted to abso- 
lute magnitudes (using the distance modulus of 30.83, Tonry 
et al. 2001 1, the globular clusters in our sample span the range 
-11.6 < M R < -8.1 in Rand -11.25 < My < -7.6 in V (for 
V- R — 0.47, Dirsch et al. 2005 1. Assuming a typical mass- 
to-light ratio of M/Ly - 2, one finds that the masses of the 
brightest clusters are of the order 7 - 8 x 10 6 Mq. For compari- 
son: the absolute visual magnitude of a> Centauri (NGC 5139), 
isM v = -10.29(Harris 1996) and its mass is about 5x10 6 Mo- 



Fig. 5. Spatial distribution of the globular clusters overlaid on 
a 33' x 33' DSS-image centered on NGC 4636. North is to the 
top, East to the left. GCs with velocities larger than the sys- 
temic velocity of NGC 4636 are marked by crosses. Circles de- 
note GCs with negative relative velocities. The concentric cir- 
cles have radii of 2, 7, and 15 arcminutes) 



5.3. Spatial Distribution 

As can be seen in Fig. [5] our target clusters are not uniformly 
distributed with respect to the center of NGC 4636. At the time 
of the mask preparation - before the photometric data set was 
fully analyzed - the projected globular cluster number density 
was overestimated. Some fields were placed at larger projected 
distances (see Fig. ^ - at the expense of a more uniform az- 
imuthal coverage. Yet in total, only 15 of the 174 GCs in our 
sample are found at projected radii larger than 7:5. Thus, any 
measurement of the line-of-sight velocity dispersion beyond 
this radius is clearly doubtful. Note that the radius at which 
the power-law index of the projected globular cluster density 
drops (cf. section fl~2l lies in the range 7'to 9'. 
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Fig. 6. Top: Color magnitude diagram for the globular clus- 
ters. Plotted are C-R vs. R. Bottom: Color distribution for the 
GCs (bin width = 0.05 mag). In both panels, the dotted line at 
C-R = 1.55 indicates the color dividing blue from red clus- 
ters (Dirsch et al. 20051 . 
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Fig. 7. Globular cluster luminosity distribution (bin width = 
0.2 mag). The dotted line at 23.33 indicates the turn-over mag- 
nitude of the GCS. Histograms for blue (solid grey) and red 
clusters. 

5.4. Velocity Distributions 

In the left panel of Fig.[8]the velocities are plotted against the 
R magnitudes. As expected, they are found to be concentrated 
towards the systemic velocity. The brightest clusters however, 
tend to avoid the systemic velocity. A larger sample is needed 
to confirm whether this effect is real or due to low-number 
statistics. The right panel of Fig.|S]shows the velocities plotted 
against the C-R colors of the clusters. Again, no clear trend is 
visible, although one notes that the five clusters with the high- 
est velocities are blue. The top right panel of Fig.|9]plots all 
GC radial velocities versus the projected radii in arcminutes. 
This figure once more shows how few velocities were obtained 
for large radii. The more important observation, however, is 
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Fig. 8. Left: Radial velocity as a function of apparent R- 
magnitude. Right: Radial velocity of the GCs plotted vs. the 
C-R -color. 



that the scatter of the velocities about the systemic velocity is 
rather small at a radial distance of about 2'-4'. As will be shown 
later (see section 0, this translates directly into a decrease of 
the projected velocity dispersion in this radial interval. On the 
other hand, only relatively few GCs fall into this radial range. 
The middle and lower right panels of Fig.|5]show the same dia- 
gram for blue and red clusters, respectively. In the radial range 
around 2'-4', the scatter about the systemic velocity is larger 
for the blue than for the red clusters, hence the radial behavior 
of the velocity dispersions for the subsamples will turn out to 
be rather different (see section^). 

6. Rotation 

Due to the inhomogeneous azimuthal coverage of our sample it 
is obvious that statements concerning rotation have to be con- 
sidered with caution. The diagnostic diagram which we have 
to analyze is a plot of radial velocities vs. the position angle. 
Cote et al. (2001 1 give a useful discussion of the relation be- 
tween the intrinsic and projected rotational velocity field of a 
spherical system and we do not repeat that here. If the intrinsic 
rotation velocity field is stratified on spheres, and the galaxy is 
not seen pole on, we measure radial velocities that depend sinu- 
soidally on the azimuth angle. Therefore, we fit the following 
relation : 



VrQS) 



+ A sin(0 - O ), 



(1) 



where v r is the measured radial velocity at the azimuth angle 
0, Vsys is the systemic velocity, and A the rotation amplitude. 
Figure [TO] shows the radial velocities of three samples ver- 
sus the azimuth angle which goes from North past East. The 
uppermost panel is the full sample, then follow the blue and 
the red clusters. A very marginal rotation signal is present in 
the full sample with the parameters A - 28 + 18kms _1 and 
©o = 63 ± 43. The blue sample does not show a significant 
signal (A = 11 ± 27kms~ 1 ), while the red sample shows sig- 
nificant rotation, which leaves its imprint on the whole sample. 
The signal and its significance depends on the color dividing 
blue and red clusters. The strongest signal is found for clusters 
redder than C-R = 1.6. In this case, we find for the entire red 
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Fig. 9. The left panels show the histograms of the GC velocities for the entire, the blue and the red subsample, respectively. In 
each panel, the grey histogram shows the distribution of the corresponding error-selected (Av < 35kms _1 ) sample. The dashed 
line indicates the systemic velocity of NGC 4636, and the bin width is 70kms~ 1 . The right panels show the radial velocities 
versus projected galactocentric distance for all, the blue and the red clusters. The crosses in circles denote data points with 
Av < 35kms~', i.e. they form the error selected samples. Plain crosses denote data points with larger velocity uncertainties. 
Again, the dashed line marks the systemic velocity of NGC 4636. 



sample, A - -87 ± 18 kms and ©o = 60 ± 22, which means 
a rotation around the minor axis with the southeastern region 
approaching. Both the amplitude and the position angle remain 
within the uncertainty limits for subsamples selected in radial 
bins. Selecting red clusters with projected distances below 5 ar- 
cmin, yields A = -105±36kms~' andQo = 72 + 24. Radii less 
than 4 arcmin give A = -114 ± 44kms~' and © = 72 + 25. 
However, a confirmation of the rotation of the red clusters is 
desirable given the relatively low number statistics and the in- 
homogeneity regarding radius and position angle. Attempts to 
interpret the rotation, for example, as a dynamical sign for a 
merger event signs seem therefore premature. In any case, it 
may be that rotation influences the velocity dispersion of the 
red clusters, making their use as dynamical tracers more com- 
plicated. 

7. Velocity Dispersion 

The top left panel of Fig. [9] shows the velocity distribution 
of the entire velocity sample. The bin size was chosen to be 
70 km s~ , which is larger than the mean uncertainty. The mean 



velocity of 906 ± 16 km s is in excellent agreement with the 
systemic velocity of NGC 4636 we determined from the galaxy 
spectra (cf. section lAlt . The shape of the distribution appears 
Gaussian (the Anderson-Darling test for normality yields a p- 
value of 0.81 (Stephens 1974)). This is what one would expect 
for a kinematical sample that is close to isotropic and in dynam- 
ical equilibrium (see e.g. Binney & Tremaine 1987 1. The most 
important observable which we want to extract from our data 
is the projected velocity dispersion and its dependence with ra- 
dius. Throughout this analysis, it is assumed that our velocity 
measurements are drawn from one normal distribution. To cal- 
culate the velocity dispersions presented in the following sec- 
tion, we employed the maximum-likelihood dispersion estima- 
tor presented by Pryor & Meylan ( 1993 1. We fixed the systemic 
velocity of NGC 4636 to v sys - 906 kms 1 . Then, the velocity 
dispersion <x is calculated (by iteration) according to: 
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Fig. 10. This plot shows radial velocity vs. the position angle 
(North past East). The upper panel plots the entire sample, 
the middle panel the blue clusters and the lower panel the red 
clusters. A marginal rotational signal is detected in the entire 
sample, which is due to the stronger signal of the red clusters, 
which seem to rotate around the minor axis with an amplitude 
of about 90 km s . The blue clusters do not rotate. 



where the sum is taken over all velocities and the £, are the un- 
certainties of the individual velocity measurements. The uncer- 
tainty of the resulting velocity dispersion is computed accord- 
ing to the expression given by Pryor & Meylan. Calculating 
the dispersion for all clusters yields: cr^u - 203 ± llkms -1 . 
The velocity dispersions for the blue and red cluster popula- 
tion are cr Uus , = 202 ± 15kms _1 and cr red = 199 + 17kms _1 , 
respectively, i.e. they are equal within the uncertainties. In or- 
der to study the radial dependence of the velocity dispersion, 
we chose to divide our samples radial bins of 1'.5 width. The 
results are listed in Tableland plotted in Fig. II II 

7.1. All Clusters 

The upper left panel of Fig. II II shows the radial dependence 
of the velocity dispersion for all 170 GCs within R < 9!0, 
For comparison, the line-of-sight velocity dispersion measure- 
ments of the central region from Bender et al. ( 1994 1 are shown 
as diamonds. These authors used a long-slit spectrum (with the 
slit positioned along the major axis) to derive a central velocity 
dispersion of cr c = 200 + 10 km s and for the average veloc- 
ity dispersion measured along the major axis (R < 32"), they 
found <r ma j = 200±30 km s . Within the uncertainties, our data 
points agree with the Bender et al. measurements. Albeit con- 
sistent with a constant velocity dispersion of about 200 km s , 
this plot suggests a decrease of cr(R) by about 50 km s (from 
~ 220 to ~ 170kms ) in the radial distance range of roughly 
2'- 4'. To check whether this behavior is merely an artifact 
produced by velocity measurements with large uncertainties, 
we next consider an error selected sample (Av < 35kms _1 ) 
(cf. Fig. [9] lower panel, to see which data points are omitted). 
The plot for this smaller (120 GC velocities) sample is shown in 
the upper right panel of Fig. ^2 Apparently, the decline is not 
due to large individual uncertainties. As will be shown next, 
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Table 4. Projected velocity dispersion as function of projected 
galactocentric radius. The upper section shows the values for 
all, the blue and the red clusters. The dispersions for the cor- 
responding error-selected samples are listed in the lower sec- 
tion. The data are shown in Fig. \H\ The values in brackets are 
those obtained when omitting the two (blue) GCs with veloci- 
ties above 1400 km s~'. 



this pattern is the result of the different behavior of the blue 
and red subpopulations. 



7.2. Blue and Red Clusters 

The middle and lower panels of Fig. \^\ show the radial de- 
pendence of the velocity dispersion for the blue and the red 
subpopulations, respectively. The values derived for the er- 
ror selected (Av < 35kms _1 ) samples which consist of 56 
blue and 56 red clusters are shown in the right panels. While 
the blue clusters show a constant velocity dispersion of about 
200 km s" 1 , the velocity dispersion of the red clusters drops 
from a value of ~ 240kms~ 1 to ~ 170kms~' at a radius of 
about 2'-3'. From this radius on, out to at least 5' 5, the ve- 
locity dispersion apparently remains at the constant, low value. 
Considering these results for blue and red clusters, one can at- 
tribute the drop seen for the whole error selected sample to the 
red clusters. However, the number of red clusters is rather high 
in the innermost radial bins, giving the large dispersion at small 
radii a high degree of confidence. On the other hand, the lower 
velocity dispersion at larger distances may be uncertain for the 
individual bins, but that the overall velocity dispersion is lower 
than for the innermost bin can be stated with high confidence 
as well. This is a new feature which has not been seen before 
in any of the few kinematic GC samples of other galaxies. We 
offer a simple explanation in the context of the spherical Jeans 
equation in the next section. 
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Fig. 11. The left panels show the velocity dispersion versus projected galactocentric radius for the whole, the blue and the red 
cluster sample, respectively. The right panels show the same for the corresponding error-selected Ac < 35km s~ samples. The 
velocity dispersions were calculated for radial bins of 1'5 width, and the number of velocities used for the measurement is 
shown below each data point. The values plotted here are listed in Table@] In all panels, the diamonds indicate the dispersion 
measurements from Bender et al. 1994 The open circles show the dispersion values at 6'that are obtained when omitting the two 
(blue) GCs with velocities above 1400 km s . 



8. The Mass Profile 



8.1. The Mass Profile from the Jeans Equation 

The collisionless Boltzmann equation is used to derive the 
spherical, non-rotating Jeans equation (see Binney & Tremaine 
1987 for a full treatment). 



+ 2 ((r) cr r (r) 

dr r 



-£(r) 



G ■ M{r) 



(3) 



The line-of-sight velocity dispersion cri 0S (R) of the globular 
cluster system and its radial dependence is used to estimate the 
mass profile of NGC 4636. Given that the number of spectro- 
scopically observed GCs is still far too low for an axisymmetric 
treatment, we will assume spherical symmetry. This is justi- 
fied since the deviation from sphericity is modest in the inner 
parts of the galaxy (Dirsch et al. 2005 1 where most of our clus- 
ters are found. A further argument is the high M/Lr found by 
Kronawitter et al. ( 2000 1 (see !8.3> Given that a flattening in the 
line of sight leads to an underestimation of M/L (Magorrian & 
Ballantynes 2001 1, it is hard to imagine that such a high value 
is still underestimated. 



with 



fi=l 



Ofl 



Here, r is the radial distance from the center and I is the spa- 
tial (i.e. three-dimensional) density of the GCs; cry and erg are 
the radial and azimuthal velocity dispersions, respectively. (3 is 
the anisotropy parameter, and G is the constant of gravitation. 
To derive the enclosed mass M(r), one needs to know /3. In 
other words: the potential and the anisotropy are degenerate in 
the spherical approximation. Information about the anisotropy 
must be inferred from higher moments of the velocity distri- 
bution, for which one needs much larger samples than ours. 
Merritt & Tremblay ( 1993 1 used M 87 as an example and esti- 
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mated that several hundred, perhaps thousand velocities would 
be required to detect a velocity anisotropy. 

8.2. The Jeans-Analysis: Principle 

For our analysis we use the expressions given by e.g. Mamon 
& Lokas (2005 1 (see also van der Marel ( 1994 1). Given a mass 
distribution M(r), a three-dimensional mass or number density 
of a tracer population t(r), and a constant anisotropy parameter 
B, the solution to the Jeans equation (Eq.[5} reads: 



Z(r)<r r {r) 



Gj ^)M(s)-(-j As 



(4) 



This expression is then projected using the following integral: 
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where I(R) is the projected number density of the tracer popu- 
lation. 



8.3. Luminous Matter 

To assess the mass distribution of the luminous (stellar) compo- 
nent, we proceed as follows: The photometric study by Dirsch 



et al. (I2005J provides a surface brightness profile in Kron- 
Cousins R. It reads: 



/*(*) = -2.5 log Ml + 



R 



+ a 2 1 + 



R 



(6) 



with a i =3.3xl0~ 7 ,.Ri =0.11',ai = 2.2,a 2 = 5.5xl0~ 9 ,.R 2 = 
8.5', and a 2 - 7.5. To obtain the luminosity density [Lq/pc 3 ], 
the intensity profile I{R) - 10 -0 ' '^ has to be deprojected. 
This was done using the Abel deprojection integral (Eq.0. For 
the cutoff-radius r t , we chose a value of 2000", roughly twice 
the distance of the last photometric datapoint. 
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dl(R) dR 
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Kr t ) 
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Further assumptions are spherical symmetry, and a distance of 
15Mpc. Moreover, we have to take into account that we use 
an R-surface brightness profile: in order to translate it into so- 
lar luminosities, we need the absolute magnitude of the Sun 
in this filter. We use a value of V -R - 0.370 (Bessell et 
al. I1998J for the Sun, i.e. M R © = 4.46. The deprojected lu- 
minosity density was then (numerically) integrated to yield a 
cumulative luminosity (in units of L ). The latter was fit by 
a sixth-order polynomial £ m, ■ R' (R in parsec) with the co- 
efficients: m = 6.38 x 10 8 , mi = 3.22 x 10 6 , m 2 = -70.0, 
m 3 = 5.62 x 10~ 4 , m 4 = 3.54 x 10~ 9 , m 5 = -8.92 x 10~ 14 and 
m.(, = 4.22 X 10~ 19 . The polynomial is a good representation in 
the range 2 < R < 70 kpc. 

To convert the cumulative luminosity into a mass profile, the 
stellar mass-to-light ratio is required, i.e. M/Lr. Kronawitter 
et al. (2000 1 quote a value of M/Lb = 9 for a distance of 
20.5 Mpc. Because M/L is inversely proportional to the dis- 
tance, one has to apply a factor of 1.36 to get M/L s = 12.2 at a 



distance of 15 Mpc. This value is suspiciously high. Although it 
formally fits old, metal-rich populations with steep mass func- 
tions, Bell et al. (J2003 > quote M/Lb - 10 as an upper limit 
found for their large 2MASS-sample. According to their rela- 
tion of (B-R) versus M/L B , NGC 4636 would have an M/L s of 
approximately 6. Further uncertainties are the distance (Dirsch 
et al. 2005 1 and the tangential anisotropy (Kronawitter et 
al. 12000 1 of NGC 4636. Obviously, there is some freedom in 
adopting the mass-to-light-ratio of the stellar population. First 
we assume M/Lb - 12.2, according to Kronawitter et al., but 
we shall also consider other values. In order to transform this 
mass-to-light ratio into an M/Lr value, one has to apply the 
factor of io-° A (( B - R )^-( B -^e\ With (B - R) 46 36 = 1-6 and 
(B - R)q = 0.97, a factor of 0.56 results, hence M/L R = 6.8. 
Further, we make the assumption that M/Lr of the stellar light 
is constant for the entire galaxy. This is justified since the color 
gradient of NGC 4636 is relatively shallow (A(C - T) <* -0.1 
in the range Q'.5 - 4'), and the stellar populations seem to be 
well mixed (Dirsch et al. 20051. 

For our modeling (cf. Sect. 18. 6L we assume that the stellar 
mass is contained within 70 kpc, yielding M stars = 4.2 x 10 11 
and 2.6 x 1O U M for M/L R = 6.8 and 4.3, respectively. The 
lower value shall illustrate the dependence on M/L. We adopt 
M/Ly =5.5 from Loewenstein & Mushotzky (2003 1 which 
results in M/L R - 4.3. 

8.4. Blue GCs as Tracer Population 

Since the projected number density N{R) of the blue globu- 
lar clusters follows a double power-law, the de-projection is 
somewhat tricky. We proceeded as follows: following Zhao 



(8) 



1997) we adopt a double power-law as 



N(K) = n (-) l + (-) 



-(fil-Yl)a\ 



Capitals mean values in projection. jx,Bi,ai are the slope of 
the inner power-law, the slope of the outer power-law, and the 
width of the transition region, respectively. B is a scale radius. 
This function was fitted to the number-counts, resulting in no = 
8.02 GCs x arcmnT 2 , B = 8.86', y\ = 0.33,/?! = 7.19, a\ = 
0.54. 

This function was then deprojected using Eq.0and in turn 
fitted with the same function, resulting in no = 0.38 GCs X 
arcmnT 3 , b = 8.86', n = 1.0, B x = 7.71, a x = 0.4. 

As a check, this 3D-profile was re-projected, using the 
Abel projection operator 



I(R) = 2 



r 

Jr 



e(r) r dr 
Vr 2 - R 2 



(9) 



and found to be in good agreement with our 2D data (see 

Rg. 0. 

For simplicity, the anisotropy parameter /J is assumed to be 
constant with radius. With our limited sample we cannot put 
any constraints on possible radial variations of B. 

As a tracer population we choose the error-selected blue 
clusters. The red clusters show a more complex behavior. They 
seem to rotate, their projected number density is best described 
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Fig. 12. Projected radial number density of blue GCs. The data 
points are from Dirsch et al. (2005 1, the solid line is the a 
double power-law (cf. Eq. [S} fit to the data and the dashed 
line is the re-projected 3D number-density. It deviates only 
marginally from the 2D-fit. For comparison we also show the 



projected light profile according to Dirsch et al. (2005 1. 



by three power-laws, and the projected velocity dispersion also 
shows more complex features. Therefore, we postpone any in- 
clusion of the red clusters in the dynamical analysis until more 
data becomes available. 

But also the blue clusters pose problems. As can be seen in 
Fig. 1 1 and Fig. 9, the velocity dispersion value of the outermost 
bin is strongly affected by two outliers. Including those objects 
one obtains 227 ± 46kms~ 1 instead of 121 ± 27kms -1 . The 
high value would mean an almost constant dispersion out to 7 
arcmin. As seen later on this would be inconsistent with a sharp 
drop of the cluster system at this radius. This unsatisfactory sit- 
uation can as well only be remedied with future data. For now, 
we calculate our models for both values. This argument cannot 
be transferred easily to the red clusters. If they rotate, a constant 
velocity dispersion is plausible since we average azimuthally. 

8.5. Dark matter 



For the dark component we adopt a density profile of the form: 

£?dark = £>dark,0 1 11 + . (10) 

V '"dark / V /"dark / 

i.e. an NFW profile (Navarro et al. H9971 Bullock et al. l2"U0TV 
For this density profile, the cumulative mass M(r) reads 
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The task is now to find values for the characteristic density 
Pdark.o and the characteristic radius rd ar k resulting in a dark com- 
ponent which, when added to the luminous matter can repro- 
duce the observed projected velocity dispersion in the radial 
range of the globular clusters. 

8.6. Results 

Having chosen a value for the stellar M/L and /3, we calcu- 
late according to Eq.0]the product of radial velocity dispersion 



squared and the 3-dimensional number density. This we project 
according to Eq. [5] For each pair of M/L and /?, we vary the 
dark matter parameters {r^y, £dark,o)> unt il me projected model 
dispersion shows minimal residuals with regard to the observed 
dispersion. For simplicity, the four data points were given equal 
weights and the residuals were calculated using 



E(o"model - 0~los) 



(12) 



We defined a set of radii (r dark e {5, 10, 40, 80, 120, 240} kpc) for 
which we then adjusted the density. Table|5] lists the model 
parameters together with derived quantities such as the virial 
mass and radius, the M/Lg-values, and the dark matter frac- 
tions within 30 kpc and 7.5 kpc (the effective radius). We de- 
termined the virial radius according to the Bullock et al. def- 
inition (i.e. at i? v j r , the mean density of the halo is 337 times 
the mean mass density of the universe, po • ^m)- Note that the 
virial quantities were calculated for the dark component alone, 
with the aim of comparing their properties to those of simulated 
N-body halos (cf. Sect. 19. 1> . Thus, having omitted the stellar 
component, the virial masses are lower limits. 

Our models are displayed in Figs. [U] and [H] which show 
our analysis for the two different stellar M/L^-values of 6.8 
and 4.3, respectively. Table|5]lists the halo parameters. The left 
and right panels refer to the difference in the last data bin as 
seen in the topmost panels which show the comparison between 
our models and the data. The observed velocity dispersion to- 
gether with their statistical uncertainties are indicated. We em- 
phasize the importance of the uncertain outermost bin which 
is clearly dominated by systematics as mentioned above. The 
various anisotropic models are distinguished by different line 
styles which can be identified in the second panel which shows 
the associated dark matter circular velocities. The degeneracy 
of potential and anisotropy is nicely visible. Only the mod- 
els without dark matter, which are shown separately in Fig. 1151 
seem to be excluded. 

The third panels show the total circular velocities, i.e. the 
sum of dark and luminous matter. The stellar component is in- 
dicated with a dotted line. The circular velocity values derived 
by Kronawitter et al. are shown as circles. Since our stellar 
mass-to-light ratio of 6.8 was adopted from the Kronawitter 
et al. analysis, it is not surprising that they do not fit well in the 
caseofM/L=4.3. 

The bottom panels show the radial dependence of the cor- 
responding M/Ls-values (cf. Table[5]for a listing of mass-to- 
light ratios at R - 30 and 7.5 kpc). 

Of course, since the stellar M/L-value enters as a free pa- 
rameter, the dark matter fraction would rise with declining 
M/L. However, M/Lr = 4.3 is already low for an old, metal- 
rich stellar population. Further constraints on the M/L ratio 
have to come from a more precise analysis of the dynamics 
of the stellar body. 

In all four model families, the /? = 1 (i.e. totally radial) 
models perform worst. Besides being a poor fit to the data, 
these models require very massive ( M v j r > 1O 13 M0) and ex- 
tended dark matter halos. 

The observed GC velocity dispersion (calculated omitting 
the two extreme velocities) as shown in the left panels of 
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Table 5. Parameters of the dark matter halos. The first columns labels the models. The mass-to-light ratio M/L R and the dark 
matter fraction are calculated for two radii (7.5 kpc a* R e ff and 30kpc is the radius within which the GCs used in the analysis 
are found). From the NFW parameters of our model halos, we calculate the virial radius, mass and the concentration parameter 
(c v j r = Ryir/^dark) using the definitions of Bullock et al. (2001 1. Then we give the total mass within 30 kpc, mass-to-light ratios 
(R-band) and the dark matter fraction within 30 kpc and 7.5 kpc (effective radius). Note that the virial quantities were calculated 
for the dark halos without adding the stellar component. Thus, the virial masses presented below are lower limits. 
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Figs. ^] and [H] can be reproduced for a range of anisotropy 
parameters (-2 < /3 < 0.5). The tangential models result in 
less massive and more concentrated dark halos. The isotropic 
model has a virial radius of about 360-400 kpc (depending on 
the stellar M/L value) which is compatible to the extent of the 
X-ray emission. 

Although the viral masses of the dark halos differ by more 
than one order of magnitude, the total mass within the radial 
range probed by the GC dynamics changes by only a factor of 
1.5. 

Hence, we can put constraints on the amount of dark matter 
in NGC 4636 within 30 kpc (^ 4R eS ), while the shape of the 
halo remains largely unknown. 

For completeness we also included in our analysis the val- 
ues of <t\ s{R) obtained when including extreme velocities of 
two GCs in the last radial bin. The resulting almost constant 
line-of-sight velocity dispersion of the GCs as shown in the 
right panels of Figs.ll3landll4lcan only be reproduced with a 
very extended dark matter component. The best agreement in 
the sense of lowest residuals is found for the fi - -2 tangential 
models. The radial/? = 0.5 models provide rather poor fits, and 
the case of completely radial orbits ft - 1 (not shown) appears 
even worse. 



9. Discussion 

9.1. Comparison to Cosmological Simulations 

How do our model dark matter halos compare to cosmolog- 
ical N-body simulations? Is such a comparison meaningful? 
Particularly the less massive halos are expected to have experi- 
enced contraction during the galaxy formation process. On the 
theoretical side, the universality of dark matter halos is still un- 
der discussion (see Reed et al. 2004 and references therein) so 
that a possible disagreement may not provide a solid basis for 
deeper conclusions. We therefore refrain from lengthy consid- 
erations and only point out basic tendencies. 
For comparison with simulations we choose the study of 
Bullock et al. (2001 1. They analyzed a sample of ~ 5000 simu- 
lated halos in the range 10 1 ' - 10 14 Mq. They found (like others) 
that the concentration parameter c wlI decreases with growing 
mass, but the scatter in the halo profiles for a given mass is still 
quite large. A comparison with Fig. 4 of Bullock et al. shows 
that for our model family shown in the left column of Fig. 1131 
(i.e. M/Lr = 6.8, omitting the two extreme blue clusters), the 
/3 — halo (model 3) fits best to the simulations. The halos 
corresponding to the models with tangential bias have higher 
concentration parameters than halos of the same mass shown 
by Bullock et al., but may be still consistent. The models with 
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Fig. 13. Mass models for a stellar M/L-ratio of 6.8. The left and right panels are distinguished by the value of the last bin. 
The topmost panels show the measured velocity dispersions (cf. Fig.^2 an d Table[4} together with the modelled line-of-sight 
velocity dispersions. The second panels show the circular velocity curves of the dark matter components. Below that, the circular 
velocity curves for the sum of luminous and dark matter are shown. The circular velocity of the stars alone is shown as dotted 
line. Circles show the values derived by Kronawitter et al. 



(2000l. The mass to light-ratio as a function of radius is shown in the 
bottom panels. Different line types indicate results for different anisotropy parameters as labeled in the second panel on the left. 



radial anisotropics have small concentration parameters, but are 
as well not excluded. This is true even for the completely ra- 
dial models (which give the worst fits to the observations). The 
question therefore is whether, given the large range of permit- 
ted halo shapes indicated by the simulations, one can meaning- 
fully compare at all theoretical simulations with observations 
of our kind, not to mention the complications if dissipative ef- 
fects of galaxy formation are taken into account. All we can 
say is that our model halos do not contradict the simulations of 
Bullock et al. 



9.2. Mass Profile and Comparison with X-ray Mass 
Profiles 

NGC 4636 has recently been a target for various X-ray studies 
(e.g. O' Sullivan et al. l20031 Ohto et al. l2"0"0"3"l Xu et al. l2"UU2l 
Jones et al. 2002 1, emphasizing the highly disturbed X-ray pat- 
tern. However, mass profiles derived from X-rays are until now 
published only by Matsushita et al. (Q998l and Loewenstein & 
Mushotzky (2003 1. FigurefTol shows the various mass profiles. 



The two solid lines are our models with the lowest and the high- 
est mass within 30 kpc (models 1 and 14). The dotted line is the 
model with the highest dark matter content within one effective 
radius (model 7). The long-dashed line is the mass profile from 
Loewenstein & Mushotzky (2003 1. The total mass enclosed 
within 30 kpc is for the latter profile 1 .25 X 10 12 M , distinctly 
higher than our highest value. Furthermore, the mass profile 
is steeper than most of our models, rising as R . However, a 
very concentrated dark halo, even surpassing this X-ray pro- 
file at small radii, is not excluded. One can speculate that the 
difference at larger radii might have to do with different hydro- 
dynamical situations at small and large radii. 

The dashed-dotted line is the mass profile from Matsushita 
et al. 1998 (based on ASCA data), read off from their log-log 
diagram (their Fig. 3, solid line, scaled down to our distance of 
15 Mpc) which shows most distinctly the flattening of the mass 
profile at lOkpc radius. However, their range of mass profiles 
consistent with the X-ray analyses is rather large, apparently in- 
cluding the possibility that there is no such flattening. In their 
work, NGC 4636 is extremely dark matter dominated, the lu- 
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Fig. 14. Same as Fig.[0] but for a lower stellar M/L-ratio of 4.3. See Table|5]for the model parameters. 



minous mass falling short by a factor of 2-3 already at small 
radii. Matsushita et al. quote M/L = 8 for the stellar component 
but give neither the band nor the source of their photometry, 
so we cannot comment further on that. The agreement with the 
GC analysis out to lOkpc is reasonably good. Outside this ra- 
dius, the discrepancy is caused by the flattening of the profile 
for which we see no indication. A similar flattening of the mass 
profile has been observed with ASCA in the cases of NGC 1399 
(Ikebe et al. l!996l and in Abell 1795 (Xu et al . 120021 . However, 
it has not been seen in the ROSAT data for NGC 1399 (Jones et 
al. 119971 . Jones et al. (2002 1 interpret the high X-ray luminos- 
ity of NGC 4636 as the transient phenomenon of an enhanced 
cooling rate. Hence it is possible that the total gravitating mass 
of this galaxy cannot be accurately determined from the X-ray 
emission. An interesting point is the large radial extension of 
the X-ray emission found by ASCA. If the dark matter frac- 
tion in NGC 4636 is low then this emission might indicate a 
dark halo not related to NGC 4636 itself. One could specu- 
late about a dark substructure in the general Virgo potential, 
formerly hosting a galaxy group which formed NGC 4636 as 
the result of the merging of several galaxies. The quite abrupt 
termination of the GC system between 7'- 9'would speak in 
favor of this scenario. Stars or GCs, respectively, near their es- 
cape energy are expected to show a drop in their radial den- 
sity profile. If NGC 4636 was embedded in a deep dark matter 



potential reaching smoothly out to large radii then one would 
expect to find GCs at large radii which are still bound to NGC 
4636 and accordingly a sharp boundary would not fit into the 
picture. However, it is interesting to note that the blue clusters 
apparently do not trace the dark halo, as has been suggested 
for other galaxies, e.g. NGC 1399 (Forte et al., 2005 1. Deeper 
insight could come from a more complete sampling of GCs or 
Planetary Nebulae at or beyond a radius of 9', which will be a 
difficult observational task. 



9.3. Modified Newtonian Dynamics 

The CDM paradigm is the most widely accepted view onto 
structure and galaxy formation. Alternative concepts to explain 
the kinematics of galaxies without resorting to dark matter are 
until today minority views. One of these alternative concepts 
is Modified Newtonian Dynamics (MOND) which one finds 
to be discussed with increasing intensity, (e.g. Milgrom 1983 
Sanders & McGaugh 2002 1. As a phenomenological recipe, the 
Newtonian acceleration gN has to be replaced by the MOND 
acceleration in the case of small accelerations g through the 
relation 



gN = n(g/ao) ■ g 



(13) 
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Fig. 15. Projected velocity dispersions for models without dark 
matter. Upper panel: M/Lr - 6.8, lower panel: the same for a 
stellar mass-to-light ratio of M/Lr - 4.3. The line styles are 
the same as in Figs.ll3landll4l i.e. the solid lines correspond to 
P = 0, the long and short dashed lines to the tangential models 
P — —j, -2, respectively. The radial models/? = 5, 1 are shown 
as dot-dashed and dot-dot-dashed lines. Only the models with 
a high stellar mass-to-light ratio in conjunction with tangential 
anisotropies are marginally consistent with the measured line- 
of-sight velocity dispersions. 
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Fig. 16. Various mass profiles. This graph shows various mass 
profiles of NGC 4636. The solid lines are the models with 
the lowest and highest mass within 30 kpc, respectively (mod- 
els 1 and 14). The dotted line is the model with the highest 
dark matter content within 30 kpc (model 7). The dashed dot- 
ted line approximates one mass profile from Matshushita et 
al. 



Il998l, showing the flattening at about lOkpc. The long- 
dashed line represents the mass profile from Loewenstein & 
Mushotzky (2003). See text for more details. 
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being a universal constant and fi(g/ao) an unspecified function 
interpolating between the Newtonian and the MOND regime. 
Bekenstein & Milgrom ( 1984i presented a Lagrangian formal- 
ism of MOND. In case of a spherical mass distribution this 
formalism reduces to the above simple formula. 
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Fig. 17. This graph shows the ratio of the MONDian acceler- 
ation to an as a function of galactocentric radius. Equality is 
reached at about 30 kpc which already is the limit of our in- 
vestigation. This means that the deep MOND regime cannot be 
reached by dynamical tracers and that the scheme of interpo- 
lating between the MONDian and the Newtonian regime dom- 
inates the analysis. A stellar mass-to-light ratio of M/Lr = 6.8 
was assumed. 



Q 

O 

S 





- 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 i_ 






- 


350 


-0.5 _ 


300 


X^— = = = = — *=^=;^_-^ 




"■-•-..., 


250 




"'•■-■•-.... 






~ 



10 



15 
R [kpc] 



20 



25 



30 



Fig. 18. MOND rotation curves for different interpolation 
schemes are shown for a stellar M/Lr ratio of 6.8 (solid lines). 
The upper solid line is derived from Bekenstein's interpola- 
tion, the middle one from Famaey & Binney's, and the lower 
one from the "standard" interpolation (see text for explana- 
tion). The dashed lines indicate the derived circular velocities 
from GCs for the given anisotropy parameters (models 3,4,5). 
The dotted line is the stellar component alone with an adopted 
M/Lr = 6.8. The Famaey & Binney interpolation apparently 
works best. 



MOND proved to have an impressive predictive power (e.g. 
Sanders & McGaugh 2002 1. Many arguments have been raised 
against MOND, including its apparent arbitrariness, but until 
now no convincing case has been found which would disprove 
MOND. Meanwhile, MOND also has a theoretical foundation 



in the work of Bekenstein (2004i. Perhaps the most serious ob- 
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jection is that MOND does not remove the need for dark matter 
on cluster scales (Sanders 2003 , Pointecouteau & Silk 2005 1. 
However, ignoring pro's and contra's for the moment, we want 
to see MOND's performance in the case of NGC 4636. 

MOND works successfully in the cases of spiral galaxies 
and dwarf galaxies. Giant ellipticals are more difficult test ob- 
jects due to the fact that MOND effects start to become vis- 
ible at radii where the low surface brightness is an obstacle 
for reliable investigations of the galaxy's dynamical behavior. 



However, Gerhard et al. (2001 1, in their sample of 21 elliptical 
galaxies, found the onset of increasing M/L-values at accel- 
erations, which are a factor of 10 higher than in spiral galax- 
ies. This would remove the universal character of qq. Milgrom 



& Sanders (2003 1 counterargued that there were discrepancies 
between Gerhard et al.'s analysis and that of Romanowsky et 
al. ( 2003 ), indicating that the last word has not been spoken yet. 
Therefore, it is interesting to have a look at the MOND rotation 
curve of NGC 4636, expected from the luminosity density pro- 
file and a certain M/L-ratio. 

FigurefTTlshows the ratio of the MOND acceleration to «o 
in dependence on radial distance (for the "standard" interpola- 
tion function mentioned hereafter). An M/L R -ratio of 6.8 has 
been adopted. The graph shows that giant elliptical galaxies 
like NGC 4636 are nowhere in the deep MOND regime. Thus it 
is clear that the interpolation function play s a dominant role and 
the dynamics at small radii do not tell anything about MOND 
without knowledge of the transition behavior at high accelera- 
tions (for which there is no theory). A formula which has been 
widely used in the case of spiral galaxies is //(x) = 

(Begeman et al. H991i . 



Vo+* 2 ) 



However, Famaey & Binney (2005) find that the above in- 
terpolation does not fit the Galactic rotation curve well and pro- 
pose yU = x/(l —x) (but see Famaey & Binney for more detailed 
remarks on the shortcomings). This formula gives a smoother 
transition to the MOND regime, i.e. MOND effects show up 
earlier. These two formulas have no physical motivation other 
than giving satisfactory results. 

Bekenstein (2004), in his relativistic TeVeS theory of grav- 
itation, which includes the MOND phenomenology in case of 
small accelerations, develops for spherical symmetry the inter- 
polation /j(x) - y +4x ~ l m This interpolation is valid only for 
small and intermediate accelerations and may not be valid for 
the full range of acceleration with which we are dealing with. 

The results are shown in Fig. ED We adopted again M/Lr = 
6.8 (an M/Lr of 4.3 would in any case not be compatible with 
MOND). The dashed lines are the circular velocities derived 
from GCs with their anisotropies indicated. The dotted line 
is the Newtonian circular velocity due to the luminous com- 
ponent alone. The solid lines represent the MOND circular 
velocities with the uppermost one being the interpolation of 
Bekenstein, the middle one Famaey & Binney's, and the lower 
one the "standard" interpolation. Apparently, MOND repro- 
duces the circular velocities derived from GC kinematics quite 
well within the uncertainties. The minimal conclusion is that 
NGC 4636 does not provide a counterargument for MOND. 
The shape of Bekenstein's interpolation at small radii is prob- 
ably irrelevant, but also at large radii the deviation from our 



mass models is quite large and becomes larger with declining 
M/L. The interpolation of Famaey & Binney works seemingly 
better than the "standard" formula, but the existent data are, 
however, not sufficient to perform a critical test. A case against 
MOND could perhaps emerge (it only can be falsified, not ver- 
ified) if the suspected drop of the velocity dispersion was to 
be confirmed in combination with a clear tangential bias. The 
necessary data set must then allow also the measurement of 
higher moments of the velocity dispersion. NGC 4636 has al- 
ready a rich cluster system, is more or less spherical, and has a 
relatively isolated location. Thus one may conclude that criti- 
cal tests are difficult to perform with elliptical galaxies. But it 
is worth mentioning that so far, one of the most interesting ar- 
guments comes from X-rays. Buote et al. (2002) note that the 
ellipticities of the X-ray isophotes around NGC 720 at larger 
radii are larger than those expected from the equipotential sur- 
faces calculated from the visible mass distributions. This holds 
for both the Newtonian and the MOND case. 

10. Conclusions 

This is the first dynamical study of the globular cluster sys- 
tem of the giant elliptical galaxy NGC 4636. Several hundred 
medium resolution spectra were acquired at the VLT with 
FORS2/MXU. We obtained velocities for 174 globular clus- 
ters in the radial range 0'90 < R < 15.' 5, or 0.5 - 9 R e in 
units of effective radius. Assuming a distance of 15 Mpc, this 
corresponds to projected galactocentric distances in the range 
4 to 70kpc. Due to the sharp decrease in the GC radial num- 
ber density at 7'- 9'about 90 per cent of our clusters are found 
in the range V < R < 7'5, i.e within 4.5 effective radii. This 
implies that the dynamics beyond this radius cannot be deter- 
mined with the current data set. Within this radius, however, we 
find a roughly constant velocity dispersion for the blue clusters. 
The red clusters are found to have a distinctly different behav- 
ior: at a radius of about 3', the velocity dispersion drops by 
~ 50kms _1 to about 170kms~' which then remains constant 
out to a radius of 7'. The cause might be the steepening of the 
number density profile at ~ 3' observed for the red clusters. 
Following this line of thought, one would expect a dramatic 
drop in the velocity dispersion for both red and blue clusters at 
7'-9', corresponding to the observed change in the power-law 
exponent of the GC number density. The observational diffi- 
culty will be to pick up sufficient clusters in this region of low 
surface density. In order to derive a mass profile, we perform a 
spherical Jeans-analysis. For a given stellar mass-to-light ratio, 
a model dark matter halo, and an orbital anisotropy, we cal- 
culate the projected velocity dispersion and choose those ha- 
los which minimize the residuals to the observational data. An 
important obstacle is the uncertain outermost data bin, which 
permits a large variety of halos to fit. Generally, slightly tan- 
gential models fit better than radial ones. Adopting the high 
stellar M/L value of Kronawitter et al. (2000) (M/L s =6.8), we 
find within 30 kpc for /3 = -0.5 an M/L^-value of 12, corre- 
sponding to a dark matter fraction of 27%. But also values as 
high as 26 cannot be ruled out with certainty, depending on the 
treatment of the last bin. When a lower stellar M/L value is used 
(Loewenstein et al. 2003 1 (M/L«=4.3), the amount of dark mat- 



18 



Y. Schuberth et al.: Dynamics of the NGC 4636 Globular Cluster System 



ter that is needed to explain the observed circular velocities in- 
creases correspondingly. Massive dark matter halos (up to 76% 
dark matter can be required. Neither of the models can be ruled 
out. An important source of uncertainty is the sparse data at the 
outer radius. However, for most of the reasonable models, we 
found either less dark matter than X-ray studies with ROS AT at 
small radii or the halos are so concentrated that the total mass 
is distinctly smaller. One should consider that the claim for the 
very high dark matter fraction is based upon X-ray luminosi- 
ties that arise in a region suspected to be out of hydrostatic 
equilibrium, so one must proceed with caution. The flattening 
of the ASCA mass profile cannot be seen. It seems that the 
dark matter problem in NGC 4636 also depends to a large part 
on the uncertain M/L value of the luminous matter. A dynam- 
ical re-investigation of the inner stellar region would therefore 
be worthwhile. Furthermore, it would be of great interest to 
obtain more spectra for globular clusters in the radial range be- 
yond 7' to see whether the suspected decline in the velocity 
dispersion can be confirmed. To answer the question whether 
the NGC 4636 globular cluster systems rotates, it is desirable to 
obtain a more uniform angular coverage. We also investigated 
whether the dynamics are consistent with Modified Newtonian 
Dynamics and answer that positively. The difficulty is that the 
deep MOND regime is beyond our radial coverage. Because 
MOND does not predict the behavior of the circular velocity in 
the transition region one has to apply physically unmotivated 
interpolation formulae. The interpolation given by Famaey & 
Binney ( 2005 ) seems to work best. 

Acknowledgements. We thank the referee, Michael Loewenstein, for 
constructive criticism which led to an improved dynamical analysis. 
Y.S. thanks the DAAD for supporting her stay in Concepcion and ac- 
knowledges support from the Graduiertenkolleg "Galaxy Groups as 
Laboratories for Baryonic and Dark Matter" and a German Science 
Foundation Grant (DFG-Projekt HI-855/2). T.R. and B.D. acknowl- 
edge support from the Chilean Center for Astrophysics, FONDAP 
No.15010003. 



References 

Ashman, K. M. & Zepf, S. E. 1998, Globular cluster systems, 

Cambridge astrophysics series 30, Cambridge University Press, 

Cambridge 
Begeman, K. G., Broeils, A. H. Sanders, R. H. 1991, MNRAS, 249, 

523 
Bekenstein, J. D., Milgrom, M. 1984 ApJ, 286, 7 
Bekenstein, J. D., 2004 Phys. Rev. D, vol. 70, Issue 8, id. 083509 
Bell, E. F, Mcintosh, D. H., Katz, N., & Weinberg, M. D. 2003, ApIS, 

149, 289 
Bender, R., Saglia, R. P., & Gerhard, O. E. 1994, MNRAS, 269, 785 
Bessell, M. S., Castelli, F, & Plez, B. 1998, A&A, 333, 231 
Binggeli, B., Sandage, A., & Tammann, G. A. 1985, AI, 90, 1681 
Binney, I. & Tremaine, S. 1987, Galactic Dynamics, Princeton 

University Press, Princeton 
Bottinelli, L. & Gouguenheim, L. 1978, A&A, 64, L3 
Bullock I. S., Kolatt T. S., Sigad Y. et al. 2001, MNRAS, 321, 559 
Buote, D. A., leltema, T. E., Canizares, C. R. Garmire, G. P. 2002, 

Apl, 577, 183 
Burstein, D. & Heiles, C. 1982, AI, 87, 1165 
Cote, P., McLaughlin, D. E., Hanes, D. A. et al. 2001, Apl, 559, 828 



Cote, P., McLaughlin, D. E., Cohen, I. G., & Blakeslee, I. P. 2003, 

Apl, 591, 850 
de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H. G, et al. 1991, 

Third Reference Catalogue of Bright Galaxies, Springer, New 

York 
Dirsch, B., Richtler, T, Geisler, D. ,et al. 2003, AI, 125, 1908 
Dirsch, B., Richtler, T, Geisler, D., et al. 2004, AI, 127, 21 14 
Dirsch, B., Schuberth, Y, Richtler, T. 2005, A&A, 433, 43 
Fall, S. M. & Zhang, Q. 2001, Apl, 561, 751 
Famaey, B. and Binney, I., 2005 MNRAS, 363, 603 
Forman, W., lones, C, & Tucker, W. 1985, Apl, 293, 102 
Forte, I. C, Faifer, F & Geisler, D. 2005, MNRAS, 357, 56 
Gerhard, O., Kronawitter, A., Saglia, R. P., & Bender, R. 2001, AI, 

121, 1936 
Giclas, H. L. 1939, PASP, 51, 166 
Hanes, D. A. 1977, MmRAS, 84, 45 
Harris, W. E. 1996, AI, 112, 1487 
Harris, W.E., van den Bergh, S. 1981, AI, 86, 1627 
Idiart, T P., Michard, R., & de Freitas Pacheco, I. A. 2003, A&A, 398, 

949 
Ikebe Y, Ezawa H, Fukazawa Y et al. 1997, Nature 379, 427 
lones, C, Forman, W., Vikhlinin, A., et al. 2002, Apl, 567, LI 15 
lones, C, Stern, C, Forman, W., et al. 1997, Apl, 482, 143 
Kissler, M„ Richtler, T, Held, E. V., et al. 1994, A&A, 287, 463 
Knapp, G. R., Faber, S. M., & Gallagher, I. S. 1978, AI, 83, 1 1 
Krishna Kumar, C. & Thonnard, N. 1983, AI, 88, 260 
Kronawitter, A., Saglia, R. P., Gerhard, O., & Bender, R. 2000, A&AS, 

144, 53 
Kundu, A. Whitmore, B. C. 2001, AI, 121, 2950 
Larsen, S. S., Brodie, I. P., Huchra, I. P. et al. 2001, AI, 121, 2974 
Loewenstein, M., Mushotzky, F 2003, Nuclear Physics B Proc. 

Suppl., 124, 91 
Loewenstein, M., Mushotzky, R. F, Angelini, L., et al. 2001, Apl, 

555, L21 
Magorrian I., Ballantyne D. 2001, MNRAS, 322, 702 
Mamon, G. A., Lokas, E. L. 2005, MNRAS, 363, 705 
Matsushita, K., Makishima, K, Ikebe, Y, et al. 1998, Apl, 499, L13 
Merritt, D. & Tremblay, B. 1993, AI, 106, 2229 
Milgrom, M. 1983, Apl, 270, 365 
Milgrom, M. 1986, Apl, 302, 617 
Milgrom, M. & Sanders, R. H. 2003, Apl, 599, L25 
Navarro I.F, Frenk C.S., White S.D.M. 1997, Apl, 490, 493 
Ohto A., Kawano N., Fukazawa Y 2003, PASI, 55, 819 
O'Sullivan, E„ Vrtilek, I. M. & Kempner, I. C, 2005 ,ApI, 624, L77 
Paolillo M., Fabbiano G, Peres G., Kim D.-W. 2002, Apl, 565, 883 
Pointecouteau, E. & Silk, I., 2005, MNRAS, 364, 654 
Pryor, C. & Meylan, G. 1993, in "The Structure and Dynamics of 

Globular Clusters", ASP Conf. Ser. Vol. 50, eds. S. G. Djorgovski 

and G. Meylan, p. 357 
Reed, D., Govemato, F, Verde, L. et al. 2005, MNRAS 357, 82 
Richtler, T, Dirsch, B., Gebhardt, K, et al. 2004, AI, 127, 2094 
Romanowsky, A. I., Douglas, N. G, Arnaboldi, M. et al. 2003, 

Science, 301, 1696 
Sanders, R. H. & McGaugh, S. S. 2002, ARA&A, 40, 263 
Sanders, R. H. 2003, MNRAS, 342, 901 
Stephens, M. A. 1974 lournal of the American Statistical Association, 

69, 730 
Temi, P., Mathews, W. G., Brighenti, F, & Bregman, I. D. 2003, Apl, 

585, L121 
Tonry, I. & Davis, M. 1979, AI, 84, 1511 

Tonry, I. L., Dressier, A., Blakeslee, I. P., et al. 2001, Apl, 546, 681 
van der Marel, R. P., 1994, MNRAS, 270, 271 
Wegner, G., Bernardi, M., Willmer, C. N. A. et al. 2004, AI, 126, 2268 



Y. Schuberth et al.: Dynamics of the NGC 4636 Globular Cluster System 19 



Xu, H„ Kahn, S. M, Peterson, J. R., et al. 2002, ApJ, 579, 266 
Yoshii, Y, Tsujimoto, T., & Nomoto, K. 1996, ApJ, 462, 266 
Yungelson, L., Livio, M., Tutukov, A., & Kenyon, S. J. 1995, ApJ, 

447, 656 
Zhao, H., 1997, MNRAS, 287, 525 
Zwicky, F. 1939, PASP, 51, 36 

The appendix lists the heliocentric velocities of globular 
clusters and foreground stars. The object identification follows 
the scheme #field.#mask#slit. The coordinates are for equinox 
2000. The colors and magnitudes are from Dirsch et al. 2005. 
It is available in electronic form. 
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